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Abstract 

We recently showed that the S&P500 stock market index is well described by Tsallis non- 
extensive statistics and nonlinear Fokker-Planck time evolution. We argued that these results 
should be applicable to a broad range of markets and exchanges where anomalous diffusion and 
'heavy' tails of the distribution are present. In the present work we examine how the Black-Scholes 
derivative pricing formula is modified when the underlying security obeys non-extensive statis- 
tics and Fokker-Planck time evolution. We answer this by recourse to the underlying microscopic 
Ito-Langevin stochastic differential equation of the non-extensive process. 
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For certain stochastic systems there is an interesting connection between statistics and 
dynamics. A family of nonlinear Fokker-Planck time-evolution equations turns out to 
be solved by probability distributions which are characterized by Tsallis non-extensive 
statistics]!], 0, |3|]. Recently we exploited this connection to analyze the dynamics of the 
S&P500 stock index showing that price-change distributions had both non-extensive 
form and Fokker-Planck time evolution. We argued that the results should be appli- 
cable to the broad range of markets and exchanges characterized by anomalous (super) 
diffusion and 'heavy' distribution tails|| ||, [1- I* 1 this paper we now investigate how 
the Black-Scholes derivative pricing formula |§ is modified when the underlying security 
is described by non-extensive statistics. This is based on an analysis of the microscopic 
Ito-Langevin stochastic differential equation underlying the macroscopic nonlinear Fokker- 
Planck equation [Q, ^, [H], |ll| . 



I. NON-EXTENSIVE STATISTICS AND TIME EVOLUTION 

We begin by summarizing the probability distribution function (PDF) P(S, t) which is 
obtained using non-extensive statistics. Denote the value of a security at a trading time r by 
price(r). In the following we will take prices and times relative to the price at some arbitrary 
fixed reference time r . Thus S(t) = price(r + 1) — price(ro) is the (relative) security value 
at a (relative) trading time t. The desired form of P(S, t) is obtained by maximizing an 
incomplete information-theoretic measure equivalent to the Tsallis entropy: 

s q = -jl—(i-Jp(s,tyds), (i) 

subject to constraints on three momentsffj, [3], ^, 11, [12]. The resulting PDF is 



P(S, t ) = ^f ) { 1 + ~l)[3- S(t)] 2 Y^ . (2) 

Here q is a time-independent parameter indicating the degree of non-extensivity or equiva- 
lently the incompleteness of the information measure. Z(t) is a normalization constant, S(t) 
is the mean, and f3{t) is related to the distribution's variance by 

1 q< 5 



a\t)= / [S - S(t)] 2 P(S,t)dS = { (b - iqmt) ^ I (3) 



oo, q > 5 
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It was shown rather unexpectedly that distributions of this non-extensive form turn out 



to solve a non-linear Fokker-Planck partial differential equation H, KL 10 



~dt~ ~ ~dS FWP&t)] + 2 dS* ■ (4) 

Here F(S) = fiS is a linear driving term dependent on the market rate of return /i. Eq. 
is solved by distributions of the Tsallis form Eq. (||) if the parameters in the latter evolve in 
time according to 

S(t) = S^e^-^ (5) 

+ 2D^\2 - q) [(3(h) Z 2 ^)] ^ [>(3-«)(t-*i) - l] (6) 
Z(t)/Z( tl ) = W)IP<h)\-i. (7) 

Here t% is an arbitrary time; e.g., it could be the shortest measured interval after the reference 
time To, so that t\ — T\ — r equals, say, one minute. 

In j|] we investigated price changes in the S&P500 index. We showed that price-change 
distributions were well-described by distributions of the non-extensive form Eq. (0) evolving 
in time according to Eq. The super-diffusion and fat tails characterizing this market 
are both a consequence of a non-extensivity parameter q greater than unity. 

The nonlinear Fokker-Planck equation is a macroscopic description of how a probability 
distribution evolves in time. It is connected to an Ito-Langevin stochastic differential equa- 
tion which describes how a particular trajectory evolves^, |TIfl . The Ito-Langevin equation 
can be written in the general form 

^ = a(S,t) + b(S,t)r ] (t). (8) 

with a the drift coefficient and b the diffusion coefficient. In the stochastic term rj(t)dt = 
dW(t) is the Wiener process @. [r)(t) is a delta-correlated {{f]{t)r]{t')) = 5(t — t')), normally- 
distributed noise with unit variance {(r](t) 2 ) = 1) and zero mean ((r](t)) = 0).] 
Eq. ([|) has a corresponding Ito-Langevin equation of the form Eq. (|8D with 



a(S, t) = F{S) = fjtS, b{S, t) = ^JDP{S, t) 1 "". (9) 
Here the Fokker-Planck equation's driving term F(S) = fiS appears as a time-independent 



linear drift coefficient a. The diffusion coefficient in our case is b = J DP(S,ty q , which 



exhibits explicitly at the level of the microscopic stochastic process the statistical dependence 
of subsequent price changes on the macroscopic PDF P(S,t). That is, the memory effect 
representing correlations in time enters here simply via the diffusion coefficient. We argued in 
H that nonlinear Fokker-Planck time evolution can be expected in any stochastic system in 
which memory effects can be approximated in this simple manner as a probability-dependent 
diffusion coefficient. 



II. DERIVATIVE PRICING FOR NON-EXTENSIVE STATISTICS 



Now we turn to the question of how the Black-Scholes derivative pricing model is to be 
modified when the underlying security has non-extensive statistics and nonlinear Fokker- 
Planck dynamics. 

We can define one form of portfolio II = —G+ S§§ @, 0, g. This is short one share of 



a derivative G and long dG/dS shares of the underlying security (stock, say) S. The change 
in the value of the portfolio in a time dt is 

BC 

dU = -dG + dS—. (10) 
uo 

(The number of stock shares dG/dS is of course constant during dt.) The change dG during 
dt is given by Ito's formula ||. This results from Taylor expanding to first order in dt and 
to second in order in price change dS, and using dW{t) 2 = dt: 

= (™ + ¥^) dt ^ ds - <»> 

Using Eq. (|S|), this becomes the stochastic differential equation obeyed by the derivative 
G(S,t). 

We seek a portfolio which instantaneously earns the same rate of return r as a short term 
risk-free security (assuming no arbitrage). Then 

dU = rUdt = r (-G + S^j dt. (12) 

Substituting Eq. (|TTD into Eq. (fLOf) and equating the result to Eq. ( fT2|) gives 

^ + rS ^! + ^ (S , 1) ^. rC(s , () . (13) 
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This pricing equation is analogous to the Black-Scholes result || , generalized for an arbitrary- 
diffusion term b(s,t) The explicit dependence on the market rate of return [i has been 
replaced by the risk-free rate r. 

There is however a difficulty hidden in Eq. (|X3|) . For our case the diffusion term b(S, t) = 



V DP x ~ q depends on the probability distribution function of the underlying stock. Hence 
Eq. ( |TB| ) depends implicitly on the market's rate of return /i. To show what difficulty this 
entails, let us begin by reviewing the Cox and Ross approach to solving Eq. (|13|). 

First define a two-point function P(S, t\S', t') which obeys an equation very similar to 



Eq. 



i)P{SJlS ' J,) ^\ /l SP{S,t\S',lf)] + ~ \b 2 (S,t)P(S,t\S',t')} , (14) 



dt dS^ y ' 1 ' n 2dS 2 

but with a boundary condition P(S,t'\S',t') = S(S — S'). [Here we continue to use b 



DP(Sjt) 1 q .\ The Cox-Ross solution is based on the fact that Eq. (14) is a forward 
Chapman- Kolmogorov equation, and as such also has a corresponding backwards form || 
14]: 

dP(S,t\S',t') cl dP(S,t\S',t') 1 2 ^ d 2 P(S,t\S',t') 

dP = ^ S dS> 2 b {S ' 1 } dS* • (15) 

As an example, let us seek a solution of Eq. ( [L3f) for a European-style call option G(S, t) on 
a non-dividend-paying stock S. Following Cox and Ross, we try the form 

G(S, t) = e- riT - t] J G{S T , T)P{S t , T\S, t) dS T . (16) 

This involves the value of G at the maturity time T. For the European call option this is 

G(S T ,T) = max(5 T - X,0), (17) 

where St is the terminal stock price and X the exercise price. (At maturity the value of the 
call option is worthless if the terminal stock price is less than the exercise price; otherwise 
the value is the price difference.) Substituting Eq. (|16~|) into Eq. flUf) , one finds that P must 
solve 

dP(S T ,T\S,t) c dP(S T ,T\S,t) l u2fa ^ d 2 P(S T ,T\S,t) 

at " os r [b > t} ds~ 2 ' 

with boundary condition P(St,T\S,T) = 5(St — S). This has the form of the backwards 
Chapman-Kolmogorov equation ( |I5| ) but with \l replaced by the risk- free rate r. 
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In all the cases considered by Cox and Ross, the diffusion term b 2 was independent of 
the underlying stock's rate of return fi. Then P was the probability distribution for a risk 
neutral world, and could readily be found. 

Here however the diffusion term b 2 depends implicitly on the underlying stock's rate of 
return /i, and consequently Eq. flUf ) cannot be solved by assuming a risk neutral world. 
Assuming risk neutrality amounts to replacing fi by the risk-free rate r everywhere. Then 
Eq. ( |18|) would be identical to Eq. (|i~5|) with fi replaced by r everywhere. This could be 
solved by a Tsallis form. Unfortunately this replacement is not justified, and we have to 
turn to an alternative approach. 

In fact, we have found two variations on the Cox- Ross approach which permit straight- 
forward solutions of the valuation equation for securities with non-extensive statistics. We 
will present both. 

The first amounts to a change of variables. Define S = S(S, t) as some function of S, t. 



Then using Ito's formula as in Eq. (11), we have 

dS = adt + ~bdW (19) 

where 

, dS l l2 d 2 S dS ~,~ . 8S 
< S ^=M + 2 b lW +a dS> b ^ = b 9S- ™ 

Options written on S can be evaluated for any b which is independent of the original security's 
market rate of return \i. Let us consider the simplest case, where b = b(t) is a function of t 
only. Then the second of Eqs. (EH) can be solved for S: 



Here we have used Eqs. (H,||). 

Now we consider an option G{S,t) and a portfolio n = — G + ^S. An analysis exactly 
like that leading to Eq. (|T3|) shows that n follows the risk-free rate of return r if 

dG ~dG l 72 d 2 G 

— +rS^ + -b 2 ^ = rG. 22 

dt OS 2 OS 2 K J 

This by construction has no dependence on /i, and hence describes a risk-free universe as in 



Cox and Ross [14]. In practical applications, S can be viewed as a derivative of the original 



non-extensive security S, and G(S,t) is then an option involving S. 



A very different route to valuing options for non-extensive securities comes from con- 
verting the market's Ito-Langevin equation into a coupled process with a constant diffusion 
coefficient. This uses an idea developed for time nonhomogeneous systems 0. Consider the 
coupled process 

dS = \jiS(t) + b(S, t)y(t)] dt = adt (23) 
dy = - iy {t)dt + ^dW{t) } (24) 

where 7 is a constant. Notice that S has no diffusion term; this will permit us to solve 
a two- variable Black-Scholes-like equation (below). First, however, we need to relate the 
coupled process to the original security S. 
Eq. ( p4"D is formally solved by 

y{t)= 1 f e-rf^VOdf. (25) 

In the limit 7 — > 00 this becomes a stationary, (^-correlated Gaussian process That is, 

y(t) — > 77(2) as 7 — > 00. (26) 



Consequently as 77 — > 00 Eq. QZ5D becomes identical to our market Ito-Langevin equation 
[Eq. (|])] and hence 5 becomes the original security S. We can then analyze Eqs. for 
finite 7 and take 7 — > 00 at the end. 

Consider an option G(S,y,t). For finite 7, Eq. (|23D has no diffusion term. Consequently 
to find dG we expand to first order in t and S, and second order in y. The result is 

(dG ,dG 7 2 <9 2 G\ , , 

dG= {-m +a ^S + TW) dt + ^ dy - (27) 

Now construct a portfolio II = — G + + One readily finds that II evolves at the 
risk-free rate r if 

dG ~dG dG 1 2 d 2 G 

— + rS-^ + ry— + -L-— = rG. 28 
at dS dy 2 dy 2 

We see that again the market rate of return \x has dropped out. Thus we have obtained a 
two-variable form of the usual Black-Scholes equation. Solutions can be find exactly as in 



Cox and Ross [Ol: 



G(S,y,t) =e~ r ^ J G(S T ,y Tl T)P(S T ,y T ,T\S,y,t)dS T dy, (29) 
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where P solves the backwards equation 

dP(S T ,y T ,T\S,y,t) ~dP dP 7 2 d 2 P 

^ = —rS — - — ry— — -. 30 

dt OS dy 2 dy 2 1 ; 

Practical application of this method of valuation would require solving with a finite value of 

7, but a value large enough so that y(t) is sufficiently close to rj(t) for time scales considered. 

We can make some connection between the two approaches described above by integrating 

out y. One can define 

G(S,t) = J dyG(S,y,t) = J dS T G{S T ,T)P{S T ,T\S,t), (31) 

where 

P(S T ,T\S,t) = * J dydy T G(S T ,y T ,T)P(S T ,y T ,T\S,y,t) (32) 

G{S T ,T) = f dy T G(S T ,y T ,T). (33) 



Eq. (|3~1) ) is formally equivalent to the solution of Eq. fl22|) in the Cox- Ross form Eq. (|T6|), in 
the limit 7 — ► 00. 

We have followed the lines of Black-Scholes and Cox-Ross to develop an options pricing 
formula for securities obeying non-extensive statistics and nonlinear Fokker-Planck time 
evolution. We showed in M that the S&P500 index is well described using this approach, 
and argued there that a description in terms of non-extensive statistics can be useful for any 
market in which the stylized facts of fat tails and anomalous diffusion are pronounced. The 
pricing formula obtained here would then be useful for options based on any such market. 
As an example we obtained pricing models for a European style call option on an underlying 
asset that pays out no dividends. 

We acknowledge support from the NSF through grant DMR99-72683. 
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